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Abstract 

The formation of new shell gaps in intermediate mass neutron-rich nuclei is investigated within the relativistic Hartree-Fock- 
Bogoliubov theory, and the role of the Lorentz pseudo-vector and tensor interactions is analyzed. Based on the Foldy-Wouthuysen 
transformation, we discuss in detail the role played by the different terms of the Lorentz pseudo-vector and tensor interactions in the 
appearing of the N - 16, 32 and 34 shell gaps. The nuclei “"^O, "^^Si and are predicted with a large shell gap and zero 

^^Ca) or almost zero ("^®Si, ^"^Ca) pairing gap, making them candidates for new magic numbers in exotic nuclei. We find from our 
analysis that the Lorentz pseudo-vector and tensor interactions induce very specific evolutions of single-particle energies, which 
could clearly sign their presence and reveal the need for relativistic approaches with exchange interactions. 
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1. Introduction 

In recent years the role of the tensor interaction in nuclei 
has gained renewed interest ( 11 ]. Essentially induced by the 
Lorentz pseudo-vector (PV) pion and the Lorentz tensor (T) 
rho meson-nucleon couplings, the tensor interaction is known 
to play an important role in the binding of light nuclei, such as 
the deuteron. At this point, we must make clear the terminology 
employed throughout this paper. The name ’’tensor interaction” 
will be used for the part of the non-relativistic nucleon-nucleon 
NN interaction which behaves as an irreducible second-rank 
tensor (2]. Such a ’’tensor interaction” can originate from a 
non-relativistic reduction of PV or T meson-nucleon couplings. 
Note, however, that the PV and T couplings contribute only to 
the exchange part of the nucleon-nucleon interaction, not to its 
direct part. Therefore, relativistic approaches such as the rel¬ 
ativistic mean field (RMF) where Fock terms are dropped will 
not lead to any tensor interaction. It was however difficult to 
assign a clear effect of the tensor interaction in medium-heavy 
and heavy nuclei. As a consequence, many-body approaches 
going from the non-relativistic Skyrme and Gogny effective in¬ 
teractions to the RMF approach, see e.g. Ref. ||3t] and references 
therein, have simply ignored the presence of the tensor interac¬ 
tion, without loss of precision with respect to global properties 
of nuclei, e.g., masses or charge radii. The situation has re¬ 
cently changed when it was shown that the nuclear tensor force 
could play a very important role in the isotcmic evolution of 
some single-particle (s.p.) states SBSHiSEIillE 
13lll4l,ll3, causing the disappearance of the usual magic num¬ 
bers and the emergence of new ones in some extreme cases of 
very neutron-rich nuclei lElIBEIIl- It is now clear that the 


modelling of exotic nuclei requires the development of more 
complete nuclear effective interactions that include the nuclear 
tensor force. 

Employing the non-relativistic approaches based on Skyrme 
interactions, the improvement of the isotopic evolution of s.p. 
energies by introducing the tensor force was found to be sys¬ 
tematically correlated to a degradation of the bindi^ ene rgy 191] . 
The tensor interaction considered in Refs. 12iil [I^lii is 
a non-relativistic contact interaction. The effects of a finite- 
range form of the tensor interaction were also explored in the 
case of the finite-range Gogny force SS. It was shown that 
the finite-range tensor force has a large impact on s.p. ener¬ 
gies along isotopic and isotonic chains but its effect on bind¬ 
ing energies was not discussed. Einally, we must mention the 
finite-range M3Y interactions which contain the tensor inter¬ 
action and are satisfactory for binding energies and s.p. spec¬ 
tra lfTllfl2 . In this approach the in-medium tensor interaction 
is not much modified and resembles to a large extent that of 
the original one in the bare NN interaction. Shell model calcu¬ 
lations with Viowk low-momentum interactions have also been 
performed, including systematic comparisons with and without 
tensor interactions, illustrating the important role of the tensor 
terms [13]. 

In all the previously mentioned studies, the tensor interac¬ 
tion at work was mostly explored in a non-relativistic frame¬ 
work. In a relativistic framework, it however contributes only 
to the Eock diagrams, where not only the n-N and Lorentz ten¬ 
sor p-N couplings play an important role but the other meson- 
nucleon couplings also carry considerable tensor force compo¬ 
nents Since in the RME approaches |21,22] the Eock dia¬ 
grams are simply dropped, the relativistic Hartree-Eock (RHE) 
approach lll2ll23ll24|] becomes the only relativistic model which 
generates a tensor force. In this work, we compare the predic 
tions based on relativistic Hartree-Eock-Bogoliubov (RHEB) 
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and relativistic Hartree-Bogoliubov (RHB) 112111221] approaches, 
which can both reproduce satisfactorily global properties such 
as binding energies and radii of finite nuclei. In the RHFB 
approach, the ;7r-pseudo-vector (;7r-PV) and p-tensor (p-T) cou¬ 
plings can be taken into account while they do not contribute 
in the RHB models. Considering the Foldy-Wouthuysen (FW) 
transformation the relation between the Lorentz type and 
rank-2 tensor forces could be analyzed, as well as their respec¬ 
tive contributions to s.p. energies and shell evolutions. We in¬ 
vestigate the shell evolution of Ca isotopes, N = 16, 32 and 34 
isotones and analyze the appearance of a large shell gap in 
and making = 16, 32 and 34 candidates for being new 

magic numbers in neutron-rich nuclei. 

This paper is organized in the following way: In Sec.|2]the 
relativistic framework is briefly recalled, as well as the non- 
relativistic reduction of the Lorentz PV and T couplings. The 
results and discussions for the selected isotopes and isotones are 
given in Sec.|3] and the summary and conclusions are presented 
in Sec.|4] 

2. Lorentz pseudo-vector and tensor interactions in a rela¬ 
tivistic Hartree-Fock model 

Within the relativistic framework, Lorentz PV Ji-N coupling 
and T p-N coupling can be introduced into the effective La- 
grangian. The corresponding interaction vertices are: 

r^''(l,2) s —■ [fnfysjyd'') (la) 

(lb) 

where m„ and M denote the mass of the n meson and of the 
nucleon, respectively, and fp are the coupling constants of 
the TT-PV and p-T meson-nucleon vertices. The other coupling 
channels, namely the isoscalar scalar cr (cr-S) and vector u) (co- 
V), the isovector vector p (p-V), and the vector photon A (A-V) 
couplings are identical to those presented in Refs. \ 2^ 

In the non-relativistic limit, using the FW transformation 
the finite-range part of the one-Tr exchange potential can be 
divided into two terms: the tensor and the central potentials, Vj 
and V^, 


V^(r) = ■T2 Si2(i + — + 

1 ^ 

= -A ■!’l ■T2 0-1- 0-2 


where 5 12 is a standard rank-2 tensor operator, 

S i2(r) = 3 ( 0-1 ■ er)(o-2 • e^) - o-i ■ 0-2. 


(2a) 
(2b) 

(3) 

The FW transformation ll^ could similarly be applied to 
the Vp potential originating from the Lorentz T p-N coupling. 
One must notice that, in the non-relativistic formalism, the cen¬ 
tral and rank-2 tensor forces generally appear independently 
from each other, while here they originate from the same in¬ 
teraction vertices. The central term plays an important 
role in determining the shell structure, which is different but as 


3. Results and discussions 

In our relativistic approaches the Lorentz PV and T cou¬ 
plings, when present, cannot be switched off without degrading 
the accuracy of the model. We instead compare various La- 
grangians having, or not, the Lorentz PV and T terms as ex¬ 
plained below. For the RHFB approach, we consider the ef¬ 
fective interactions PKAl (27] and PK03 ifioll . The former is 
the most complete relativistic model to date, taking both tt-PV 
and p-T into account, while the latter does not contain the p-T 
coupling. All meson-nucleon couplings are density-dependent 
and their values have been adjusted in previous studies. Two 
density-dependent RHB models are also used throughout this 
work: DD-ME2 (H and DD-MEd Q. The DD-ME5 differs 
from DD-ME2 by the inclusion of the d-meson, which leads 
to different proton (n) and neutron (v) Dirac masses. Since all 
these modelings do not contain originally a term describing the 
pairing effects, the Gogny force DIS Oil] is employed in the 
pairing channel. DD-ME2 and DD-MEb do not have tensor 
terms (neither of Lorentz type, or of rank-2 type). Since these 
four effective Lagrangians reproduce equally well the energies 
and radii of finite nuclei, we could compare in the following the 
impact of the tensor force on more detailed quantities such as 
s.p. e nerg ies, within a spherical model which was described in 
Refs. Him. 


important as the rank-2 tensor term 01 ll l28l] . 


3.1. Lorentz pseudo-vector and tensor effects on Ca isotopes 

The analysis of the oxygen isotopes has revealed the very 
important role of the spin-isospin properties of the nuclear in¬ 
teraction and their impact on the magicity of \ TH 28 ]. The 
origin of this strong spin-isospin interaction is directly related, 
in a one boson-exchange picture, to the contribution of n and 
tensor-p mesons. In a similar way, it is interesting to explore 
the role of these mesons along the next semi-magic Z = 20 
(resp. N - 20) isotopic (resp. isotonic) chains with respect 
to the neutron (resp. proton) shell evolution. Therefore, it is 
expected that, since there exist a large number of measured 
properties, the analysis of the theoretical predictions could shed 
light on the role of these meson fields, and more precisely, of 
the tensor force. Notice that "^°Ca is spin-saturated in the proton 
and neutron shells, and the tensor interaction that contributes 
in calcium isotopes above "^°Ca is mostly active between a few 
neutron states: mostly v2p and vl/ states. Calcium isotopes 
therefore provide an ideal isotopic chain for the theoretical and 
experimental analysis of the neutron-neutron tensor interaction 
and for studying its role in the formation and evolution of neu¬ 
tron shells in medium mass nuclei. 

We first analyze the question of the possible proton level in¬ 
version in ^*Ca. The normal level sequence being {7r2ii/2, n\d 2 ,i 2 ], 
it has been suggested that the inversion of this sequence if well 
marked, could induce the formation of a proton depletion, a 
kind of bubble structure, in the neighboring "^®Ar, reducing the 
spin-orbit (SO) splitting 1^ . The energy difference = 
S 7 r 2 .si /2 “ SttWsp along the calcium isotopic chain is shown in 
Fig.IU for the selected effective Lagrangians: PKAl, PK03, 
DD-ME2 and DD-ME^. The experimental data of the rele¬ 
vant s.p. energies in both ^°Ca and ^^Ca are also reported in 
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Figure 1 : (Color online.) Energy difference = s^2si/2 ~ ^^1^3/2 “f Ca 

isotopes. The experimental values of ^*’Ca and '**Ca are taken from Ref. [3 
(up-triangle), Ref. O (down-triangle), Ref 13 (left-triangle), and Ref I37I1 
(right-triangle), respectively. 


Refs. 1134 135 LI36ll37n . Experimentally, s.p. energies are diffi¬ 
cult to determine due to the st rength frag mentation caused by 
the particle-vibration coupling 1^, 3^ Among these effec¬ 
tive Lagrangians, it appears that only the RHFB-PKAl model 
gives a satisfactory agreement with the experimental results and 
predicts an almost perfect degeneracy between the 7r2si/2 and 
nld^l 2 states, consistent with the data (see Fig.[TJ. The other 
effective Fagrangians, which do not contain the Forentz PV and 
T couplings or only partially (such as PK03 which has the PV 
coupling), fail to reproduce the experimental points, and a level 
inversion is even predicted around ^^Ca, in contradiction with 
the data. From the comparison shown in Fig. [1] we deduce that 
the tensor-p meson-nucleon coupling, which can be treated as a 
mixture of central and tensor forces, is an important ingredient 
in order to reproduce the s.p. spectra in Ca isotopes. The impor¬ 
tance of the tensor force in this context has also been stressed 
by non-relativistic approaches 114, T^]. It is worthwhile to 


notice that the two proton states nls 1/2 and Tild-ii 2 are pseudo¬ 
spin (PS) partners ll43[ 14411 . The energy difference could 
therefore be interpreted as a measure of the PS degeneracy. In 
this case, the Forentz tensor p-field is also an important ingre¬ 
dient for discussing the occurrence of PS degeneracy in "^^Ca. 

We now turn to the discussion of the emergence of new shell 
closures in exotic calcium isotopes which have been indicated 
by recent experiments. First, the high excitation energy of ^^Ca, 
compared to those of neighboring nuclei, favors a possible new 
N - 32 shell closure 14511 . Then, a confirmation was obtained 
from recent high-precision mass measurements of several iso¬ 
topes ranging from ^'Ca up to ^'^Ca ll46ll . In neighboring tita¬ 
nium and chromium nuclei, a high excitation energy was also 
observed indicating the robustness of the N -32 magic number 
in this region of the nuclear chart 114711481] . Additionally, from 
the very recent measurement of the 2| energy in ^^Ca which 
was found to be only ~500 keV below that in ^^Ca ll49n . it is 
also possible that ^^Ca could be magic as well. If confirmed 
it would be the first example of two successive magic numbers 
(N - 32 and 34) in the same element. It should however be no¬ 
ticed that the possible V = 34 magic number is less marked in 
the neighboring nuclei such as titanium where the 2| excitation 


energy was found smaller than in calcium 


The role of the different components of the nuclear interac- 



Figure 2 : (Color online.) Upper panels: Neutron s.p. spectra of Ca isotopes, 
extracted from the RHFB calculations with PKAl (plot a) and the RHB ones 
with DD-ME 2 (plot b). Lower panels: Neutron pairing gaps determined by the 
RHFB (plot c) and RHB (plot d) calculations. The experimental data (in various 
triangles) are taken from Refs. see Fig.^for the convention of 

symbols. 


tion can be analyzed from various perspectives. The mean-field 
in non-relativistic Skyrme approaches allows to analyze the re¬ 
spective roles of the central, SO and effective mass terms 114111 : 
the tensor force in M3Y-type models was found to be important 
in reproducing the experimental trends El; the comparison be¬ 
tween the relativistic and non-relativistic approaches shows that 
N ox Z - 8, 20 are suitable for fitting the tensor parameters m, 
and finally, calcium isotopes have been used as benchmark nu 


clei to illustrate the role of the three-body (TB) force 11521153 


15411 . Concerning the last study, it is worth noticing a recent 
work showing that a properly optimized chiral two-body inter¬ 
action can also describe many aspects of nuclear structure in 
oxygen and calcium isotopes without explicitly invoking TB 
force 1551]. Ju st after the measurements of neutron-rich calcium 
isotopes i4a. 14911 . it was argued from a Skyrme Hartree-Fock 
(SHF) approach that the magicity of ^^Ca and ^^Ca might be 
reproduced selecting fitted tensor parameters Ei- Pairing cor¬ 
relations have however not been considered and the magicity 
has been inferred only from the size of the s.p. gaps. It is 
known that relatively large s.p. gaps could still be overcome 
by pairing correlations leading to fractional occupations above 
the expected s.p. shell gap. As proposed in a non-relativistic 
approach based on the Skyrme force E^ . a criterion for the 
magicity should thus also refer to the occupation number above 
the s.p. gap which, in case of magicity, should be found to be 
zero. 

Fig.|2]shows the isotopic evolution of the canonical s.p. en¬ 
ergies (upper panels) and average pairing gaps (lower panels) 
from "^°Ca to ®®Ca, going over the possible new magic num¬ 
bers at A = 32 and 34, and using RH(F)B approach with the 
selected effective Fagrangians. In Figs. |2] (a) and (b) there 
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Table 1: Energy difference A£'(/, f') = £/ - Sj' (in MeV) in Ni, Ca and Si at 
N = 32 and 34. The results correspond to RHFB-PKAl, PK03 and RHB-DD- 
ME2 Lagrangians. See text for details. 


tional suppression of the neutron pairing gap at = 32 which 
is not predicted by the other Lagrangians. This large quenching 
for = 32 is an additional hint which suggests that the N = 32 


Eorce 

^E{i, i') 

A 

Ni 

Ca 

Si 

magic number should be analyzed in the light of the Lorentz PV 

PKAl 

(v 2 pi/ 2 , v 2 p 3 / 2 ) 

32 

1.51 

2.72 

0.81 

and T interactions. A weaker, but still distinct and important. 

(vl/ 5 / 2 , v 2 pi/ 2 ) 

34 

1.04 

2.45 

4.05 

quenching is also predicted by the RHEB-PKAl at A = 34, 

PK03 

(v 2 pi/ 2 , v 2 p 3 / 2 ) 

32 

1.22 

1.69 

0.68 

which may suggest as well A = 34 to be a submagic number. 

(vl/ 5 / 2 , v 2 pi/ 2 ) 

34 

-1.72 

0.77 

2.72 

It is also interesting to notice the predictions at A = 40: the 

DD-ME2 

(v 2 pi/ 2 , v 2 p 3 / 2 ) 

32 

1.58 

1.76 

0.92 

RHB models indicate a quenching of the pairing correlations at 

(vl/ 5 / 2 , v 2 pi/ 2 ) 

34 

-1.23 

1.21 

3.18 

A = 40, while the RHEB ones predict a small decrease but not 


are shown the evolutions of the neutron s.p. energies around 
the chemical potential determined by RH(F)B calculations with 
PKAl and DD-ME2, respectively. We remind that the RHB- 
DD-ME2 model contains neither the Lorentz tensor nor rank-2 
tensor force com pon ents. In Eig. |2] the experimental values 
taken from Refs. 35, 3^ 33] are also shown. It should be 
noticed that we only take the experimental data as a reference 
due to the limit of mean-field approach. In the present calcu¬ 
lations, the correlations beyond mean field 1 38, [BlidlsSIsl, 
such as the particle-vibration couplings, are not taken into ac¬ 
count. Another fact is that among the selected s.p. states some 
states such as the v 2 / 7 i /2 in ^°Ca or vl/ 5/2 in "^^Ca are highly 
fragmented. 

The general feature reflected in Eig.|2]is that, as the number 
of neutrons increases, the s.p. energies decrease due to the en¬ 
hanced mean-field potential. This trend is smooth for the RHB- 
DD-ME2, as well as for the other non-presented models RHEB- 
PK03 and RHB-DD-MEb, while we observe abrupt changes 
in the RHEB-PKAl results when the s.p. energies cross the 
chemical potential. Notice moreover that "^°Ca is spin-saturated, 
while the spin asymmetry reaches a maximum at ^^Ca. We 
have analyzed the isotopic evolution of the two-body interac¬ 
tion matrix elements V,/- between neutron valence orbits and 
we found that the couplings are quite constant for the RHB- 
DD-ME2 while for the RHEB-PKAl they are slightly more de¬ 
pendent on the number of neutrons. In particular, it is observed 
a reduction of the coupling between v 2 pi /2 and v 2 pj ,/2 states at 
A = 32 by about 40% as compared to the values at A = 30 
and 34. This reduction is associated with the emergence of a 
shell gap at A = 32. We also observe that the coupling between 
v2pi/2 and v2p2/2 states is reduced by a factor 2 with the RHEB- 
PKAl as compared to the RHB-DD-ME2, which explains why 
the gap at A = 32 is not so distinct for the RHB-DD-ME2. In 
contrast, although both Lagrangians give almost identical cou¬ 
pling strengths between the v 2 pi /2 and vl/ 5/2 states, only the 
RHEB-PKAl shows a gap at A = 34. The gaps at N - 32 and 
34 are therefore related to different parts of the nuclear interac¬ 
tion, and we will analyze it with respect to the Lorentz PV and 
T interactions in the following. 

As discussed before, the magicity is not solely related to an 
increase of the shell gap, but also to a quenching of the pair¬ 
ing correlations at magic numbers. As seen from Pigs. |2c) and 
| 3 d), the suppression of the neutron pairing gaps at traditional 
magic numbers A = 20 and 28 is confirmed by all the models 
considered here. In the RHEB-PKAl results, there is an addi- 


Pock term and to the effect of the Lorentz PV and T interac¬ 
tions. The most advanced Lagrangians presented here (PKAl 
and PK03) therefore do not predict A = 40 as a magic number 
for “Ca. 

3.2. Evolution along A = 32,34 down to the drip line 

We now analyze the contribution of the neutron-proton in¬ 
teraction to the SO splittings of v2p states. To do so, we fix the 
neutron number as A = 32 and vary the proton number from 
“Ni down to ^^Ca, and along the isotonic line the valence pro¬ 
tons are gradually removed from the 711 / 7/2 state (which is a 
7 > state). The energy differences associated with the N - 32 
and 34 shell gaps are given in Table [1] for a set of selected iso¬ 
tones going from the stability line down to the drip line; “’®^Ni, 
and Comparing the v2p SO splittings (A = 32) 

given by different effective Lagrangians, it is interesting to no¬ 
tice that the values are almost the same for “Ni (1.4+ 0.2 MeV) 
and for ^®Si (0.8 + 0.1 MeV), while they are quite different for 
^^Ca. A more systematic calculation along the N -32 isotonic 
line is shown in Pig. [3] where Pig.|3la) shows the comparison of 
the SO splittings between the models (PKAl, PK03, DD-ME2 
and DD-ME(5) and Pig. |3b) presents the detailed contributions 
of PKAL To better analyze the isotonic evolution, we present 
the results with respect to the values calculated in ^^Ca. 

As shown in Pig. [3 a), the RHEB-PKAl results show a dis¬ 
tinct enhancement of the SO splitting from Ni to Ca (~ 1.2 MeV), 
while it is less pronounced for RHPB-PK03 (~0.5 MeV) and 
much less for RHB-DD-ME2 and RHB-DD-MEb (-0.2 MeV). 
On the other side, from S to Ca the situation is changed. To¬ 
wards Ca, the RHPB-PK03 results present a 0.8 MeV enhance¬ 
ment of the v2p SO splitting (-0.5 and 0.4 MeV respectively 
for RHB and RHEB-PKAl). Here, we truncate at S to make 
the figures more readable, but the conclusions extend beyond. 
To further analyze the origin of this behavior and the role of the 
Lorentz PV and T interactions, we have plotted in Pig.|3b) the 
contribution from the Lorentz PV and T couplings in RHEB- 
PKAl, which are decomposed into the rank-2 tensor term 
and the central term V*', as compared to the other terms (Rest) 
associated with the kinetic energy term, cr-S, w-V and p-V cou¬ 
plings. We find that the Lorentz PV and T couplings, namely 
the TT-PV and p-T ones in PKAl, play a dominant role in de¬ 
termining the enhancement of the v2p SO splitting from Ni to 
Ca, while from Ca to S, the contributions of the Lorentz PV 
and T couplings are close to zero. Such difference is due to 
the fact that the contributions of the central and rank-2 tensor 
terms, namely and , add up between Ni and Ca, while 
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Figure 3: (Color online.) (a) The SO splitting of v2p along the N = 
32 isotonic chain, calculated with the RHFB (full symbols) and RHB (open 
symbols), (b) Detailed contributions to the SO splitting from the ;r-PV and p-T 
couplings, in comparison with those from the other channels. The results are 
obtained with the RHFB-PKAl Lagrangian. 


Figure 4: (Color online.) Detailed contributions from the ;r-PV and p-T 
couplings to the energy differences (a) ~ ^vlpip and (b) 

- ^vl/ 5/2 ” ^vlp 2 j 2 along the A = 34 isotonic chain. The results cor¬ 
respond to the RHFB-PKAl Lagrangian. 


they mutually compensate between Ca and S. In addition, it is 
worth noticing the large contribution of the central term to the 
SO splitting between Ni and Ca (70% of the increase), while the 
rank-2 tensor accounts for only 30% of the increase. In conclu¬ 
sion, it is illustrated in the case of the N - 32 isotones that the 
effects of the Lorentz PV and T on the s.p. energy gaps are 
certainly not reducible to the contribution of the rank-2 tensor 
because the central term plays a more important role. 

We now turn to the N -3A gap which is appearing between 
the v2pi/2 and vl/5/2 states. Similarly to the N - 32 case, we 
first refer to Table [1] where the energy differences responsible 
for the N - 3A shell gap, = Svi/ 5/2 “ ^v 2 pi/ 2 , are given 

for Ni, Ca and Si. Only the RMFB-PKAI model presents fairly 
distinct shell gap at N = 34 for ^^Ca and all the models pre¬ 
dict much enhanced gap for '^^Si. For ®^Ni, the order of the 
states is even changed in the RHFB-PK03 and RHB-DD-ME2 
results. In order to clarify the role of the Lorentz PV and T in 
the N - 3A case, we represent in Fig. |4ja) the evolution of the 
energy difference AL^^ shifted by that in ^"^Ca, while the en¬ 
ergy difference involving the SO partner of the occupied state, 
AL^.^ = fiyi- ev 2 p 3 / 2 ’ represented in Fig. Stb). The role of 
the Lorentz PV and T is not so straightforward since the states 
are not SO partners. It is found, for PKAl, that the Lorentz PV 
and T forces present tiny contributions to the splitting between 
the v2pi/2 and vl/5/2 states. As shown in Pig.|3a) the origin of 
this weakening is due to the near cancellation between the cen¬ 
tral and the rank-2 tensor components of the Lorentz PV and T 
couplings. This cancellation occurs for all isotones considered 
here, from Ni to S. 

It is interesting to notice that the states {v2/73/2, vl/5/2}, in¬ 
volving the SO partner of the occupied state of the N -3A gap, 
are PS partners and usually named vld, cf. Lig.E] The contribu¬ 
tion to the PS splitting of Lorentz PV and T terms (;7r-PV and p- 
T) as well as that coming from the other terms (Rest) calculated 
with the RHLB-PKAl model are represented in Lig.lUb). It is 
observed that the other terms (Rest) largely contribute to the en¬ 


ergy difference AL^^ , while the Lorentz PV and T contribute 
only to about 30% of the splitting. Lor the N — 3A isotones, 
we observe that the v2p SO splitting is changing rather weakly 
from Z = 28 to 20, and therefore the opening of the V = 34 gap 
for Z = 20 can be understood, to a large extent, from the evolu¬ 
tion of the PS splitting. The shell gap at A = 34 can therefore 
be interpreted as a manifestation of a strong isospin-dependent 
PS splitting in which the Lorentz PV and T terms have only a 
weak impact. 

It is worth comparing our results to very recent experimen¬ 
tal analyses of shell gaps in ^‘’Ar, ^"^Ti 1^ . Using 

RHLB-PKAl, we observe as well that the N -32 subshell gaps 
in ^°Ar, ^^Ca, ^'^Ti are similar in magnitude, see Pig. [3ja), and 
that the A = 34 subshell closure in ^^Ar is larger than in ^"^Ca, 
see Pig. Ilja). In addition, the shell gaps continue to increase 
between ^^Ar and ^°S, see Pig. EJa), and this tendency is con¬ 
firmed for the next A = 34 nucleus which is '^^Si. RHLB-PKAl 
predicts a large gap for the drip line nucleus '^^Si (~4.0 MeV) 
and a small pairing gap, while the other Lagrangian predict a 
small gap, see Table [1] 

In conclusion of this part, the effect of the Lorentz PV and 
T couplings (tt-PV and p-T) on the N - 32 shell gap is shown 
to be dominant at variance with the A = 34 one. The origin 
of these shell gaps are indeed different: the N - 32 shell gap 
is related to the v2p SO splitting while the N — 3A shell gap 
can be related to the vld PS splitting. The prediction of two 
successive magic numbers in Ca isotopes (^“’^"^Ca) is not related 
to the same origin. 

3.3. Magicity ofN — 16 in neutron-rich isotopes 

Based on the evolution of neutron separation energies and 
cross sections in light nuclei, it was originally proposed A = 16 
to be a new magic number lying between the usual A = 8 and 
20 for N, O and P nuclei 1591] . In these nuclei, the s.p. gap at 
A = 16 is occurring between the vldi /2 and v2ii/2 states. Sev¬ 
eral experimental analyses have concluded as well that A = 16 
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Ca Ar S Si Mg Ne O Ca Ar S Si Mg Ne O 



Figure 5: (Color online.) (a) Energy difference = £vld^i 2 ~^v 2 sii-^ along the 
N = \6 isotonic chain, calculated with the RHFB (full symbols) and RHB (open 
symbols). The experimental value (triangle) of 1^ is also displayed as a 
reference, (b) Detailed contributions to from the ;r-PV and p-T couplings, 
in comparison with those from the other channels. The results are obtained with 
the RHFB-PKAl Lagrangian. 


could be a magic number. For instance, from proton knock¬ 
out reactions using beam, the s. p. g ap between these states 
was deduced to be 4.86(13) MeV ibOtl . Another conhrmation 
came also from the measurement of a relatively high excitation 
energy and a small (32 value in 116111 . These results also in¬ 
dicate that has a spherical character. 

Our predictions are shown in Fig. |5] where we represent the 
energy difference between the vld^ii and v 2 si /2 states, namely 
= evW 3/2 “ £v 2 ii/ 2 j along the A = 16 isotonic chain from 
Z — 2Q down to 8. Here, we have used as pairing interaction 
the original Gogny force DIS modihed by an adjusted strength 
factor f — 1.10, so as to reproduce the odd-even mass differ¬ 
ences of oxygen isotopes 16211 . In Fig. |3a) we compare the 
predictions of different effective Lagrangians, PKAl, PK03, 
DD-ME2 and DD-ME5, with the experimental data in 
The effective Lagrangian PKAl predicts a large s.p. gap from 
Z = 8 up to 20, as well as a s.p. gap in close to the ex¬ 
perimental value. At variance with this prediction, the three 
other effective Lagrangians (PK03, DD-ME2, and DD-ME5) 
predict a strong reduction of the s.p. gap when going from 
Z = 8 to 20, and the s.p. gap predicted for is well be¬ 
low the data. In order to analyze the role played by the Lorentz 
PV and T couplings on the s.p. gap evolution, we present in 
Pig. |5lb) the contributions from the central and rank-2 tensor 
terms of the Lorentz PV and T couplings, as well as the other 
terms (Rest) for the RHPB-PKAl model. Lor Z = 20 down to 
8 , the Lorentz PV and T couplings contribute about 30-50% of 
the s.p. gap, and it is interesting to note that the contribution 
from the central terms of the Lorentz PV and T couplings are 
dominant and the effect of the rank-2 tensor coupling is quite 
weak. The latter observation contradicts previous studies based 
on non-relativistic approaches using Skyrme and Gogny forces 
where a stronger rank-2 tensor force is necessary to reproduce 
s.p. evolutions Q. The contribution of the Lorentz PV and T 
tensor couplings (yellow region in Pig.|5]) for N - 16 is rather 


constant along the isotonic chain, and the observed Z-evolution 
is mostly related to the ’’Rest” terms and seems to be quite im¬ 
pacted by shell effects. A general feature observed with our 
relativistic models is that the shell gap is decreasing from O to 
Si, despite an increase around Z = 14. 

4. Summary and conclusions 

In summary, the formation of new shell gaps in neutron- 
rich nuclei is investigated within the RHPB theory and the role 
of the Lorentz PV and T interactions is studied in detail by com¬ 
paring different Lagrangians with or without such terms. Based 
on the most complete RHPB-PKAl effective Lagrangian, we 
conhrm and predict that "^®Si, and are the magic 

neutron-rich nuclei. In the case of ^^Ca, the role of the Lorentz 
PV and T components is determinant, while it is less impor¬ 
tant for and negligible for ^"^Ca and ^*Si. Analyzing shell 
evolutions along the N - 32, 34 and A = 16 isotonic chains, 
we observe that the global variation of the s.p. energies is due 
to the isoscalar component of the effective Lagrangian, while 
the more specihc evolution of SO and PS partners is related to 
the Lorentz PV and T couplings (tt-PV and p-T). Based on the 
PW transformation, we analyze the role played by the central 
and rank-2 tensor terms of the Lorentz PV and T couplings in 
the formation of the A = 16,32 and 34 shell gaps. It is shown 
that those terms drive the distinct enhancement of the N - 32 
gap along isotonic chains going from ®°Ni to ^^Ca, while the in¬ 
crease of the A = 34 gap from ®^Ni to ^^Ca, and "^^Si, is mainly 
due to the other channels. The shell gap at A = 34 can be also 
interpreted as a manifestation of a strong isospin-dependent PS 
splitting. Pinally, we also observe that the Lorentz PV and T 
couplings are not likely to support the appearance of A = 40 
magicity in ®°Ca. 

In this study, we have illustrated in intermediate mass nuclei 
the very important role played by the Lorentz PV and T inter¬ 
actions which cannot be simply reduced to their rank-2 irre¬ 
ducible tensor contribution. We have indeed observed a strong 
interplay between the rank-2 tensor and the central terms orig¬ 
inating from the Lorentz PV and T interactions. While these 
two terms generally appear separately in non-relativistic nu¬ 
clear models, they originate from the same interaction vertex 
in the relativistic ones. In the intermediate mass nuclei that we 
have explored, the s.p. energies are shown to be very much 
impacted by the Lorentz PV and T interactions which impose 
some very specihc behavior along isotonic chains. Even their 
weak inhuence results from the cancellation of two of their con¬ 
stituent terms. The effect of the Lorentz PV and T couplings — 
which are purely relativistic terms — in s.p. energies is there¬ 
fore very different from the rank-2 tensor force. We conclude 
from this study that, while both relativistic and non-relativistic 
models for the nuclear interaction could equally well reproduce 
global properties of hnite nuclei, such as binding energies and 
radii, the detailed evolution of s.p. energies could potentially 
sign relativistic effects such as the Lorentz PV and T interac¬ 
tions. This statement requires however further analysis based 
on other modelings of the nuclear interaction, but we propose 
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here a clear relation between Lorentz PV and T relativistic ver¬ 
tices and s.p. evolution. 
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